Secondary electron cascade in attosecond photoelectron spectroscopy from metals 
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Attosecond spectroscopy is currently restricted to photon energies around 100 eV. We show that 
under these conditions, electron-electron scatterings, as the photoelectrons leave the metal give 
rise to a tail of secondary electrons with lower energies and hence a significant background. We 
develop an analytical model based on an approximate solution to Boltzmann's transport equation, 
to account for the amount and energy distribution of these secondary electrons. Our theory is in 
good agreement with the electron spectrum found in a recent attosecond streaking experiment. To 
suppress the background and gain higher energy resolution, photon sources of higher energy could 
be advantageous. 
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The time-resolved study of metal surfaces using femto- 
and attosecond pulses is currently developing as a sepa- 
rate research topic [E 0, H, 0; IE @1- It is now possi- 
ble to perform experiments with a very high temporal 
resolution, allowing time-resolved investigations of elec- 
tron dynamics. These include collective dynamics such as 
plasmon creations 0, Q > electronic screening within the 
metal [j| and image-charge states [To! • When combining 
XUV pumps and infrared probes, it is possible to study 
in real time the Auger decay from atoms adsorbed on a 
metal surface Q and using the streaking technique fu\ , 
to resolve a delay of only 110 asec between photoelectrons 
emerging from the conduction band and core levels 0. 

A series of theory papers study the sh ort-p ulse interac- 
tions with metal surfaces 12, .Progress m 



this field is challenging since field-free, solid state meth- 
ods are not readily extended to account for the highly 
nonlinear and explicitly time-dependent dynamics intro- 
duced by the short pulses. Even when it is possible to 
regard the interaction in a one-electron model, many dif- 
ferent secondary processes may blur the direct process. 
Hence, development of models that capture the essen- 
tial physics and are still computationally tractable are 
called for. Here, we develop a simple, analytical descrip- 
tion of the electron-electron scattering in the metal and 
show that this leads to a large background of electrons 
~ E~ 2 . This energy-dependent background is a result 
of the spectral broadness of the attosecond pulse and 
leads to a reduced energy resolution for any signal from 
electrons below the conduction band. To suppress the 
background, high-energy photon sources should be used. 

Our theory is based on the fact that fast primary elec- 
trons created by the attosecond pulse, will produce sec- 
ondary electrons by collisions with electrons bound in the 
metal, and initiate a cascade as they move through the 
metal. The primaries and secondaries lose energy until 
they escape, if sufficiently energetic, or sink back into 
the Fermi sea. The secondaries are responsible for the 
low-energy tail in the spectrum. To find the spectrum, 



including the above mechanism, implies the solution of a 
full dynamic many-body problem. Instead of persuing a 
quantum solution, we aim at a semi-analytical approach, 
where the energy distibution is determined by the solu- 
tion of a transport equation for the electrons. To formu- 
late the transport equation, we use that the main way the 
electron loses energy in the energy range of interest is by 
collision with conduction electrons; the electron-phonon 
coupling is weaker 17| . In a steady-state regime, and 
based on neutron transport theory [l8( , casade processes 
pertaining to this situation were considered many years 
ago for an electron beam impinging normally to a metal 
surface 17). In the present scenario, however, the pulsed 



light sources introduce an explicit time-dependence in the 
problem, and our starting point is the time-dependent 
Boltzmann equation [l^ | 
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where $(r,fi, E,t)drdftdE is the electron distribu- 
tion function, i.e., the number of electrons in [r;r + 
dr], [ft; ft + dfl], [E; E + dE] at time t, r is the posi- 
tion within the metal, ft = v/v is a unit vector in the 
direction of the electron velocity v, and E is the energy. 
The LHS in ([TJ) , represents the time-rate of change of the 
electron distribution function in the direction ft, includ- 
ing the effect of the force F(t) from the streaking near-ir 
fs pulse. On the RHS the first term describes the loss 
of electrons due to scattering and possibly capture, and 
X(E) is the mean free path. J is the source term re- 
sulting from the attosecond pulse, giving the number of 
electrons emitted per unit time at r and t with energy E 
in direction ft. The factor u'$(r, ft', E',t)/X{E') in the 
last term describes the number of collisions per unit time 
at r occuring to electrons with energy E' and direction 
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SI', finally, S(E, SI; E, SI') is the probability that, given 
an electron at E',Sl', one will be found at E,Sl after a 
scattering. 

First, we disregard the fs, ir-pulsc and focus on the 
influence of electron propagation on the number of elec- 
trons passing per unit surface area induced by the at- 
tosecond uv pulse alone. We therefore leave out the force 
term and integrate ([1]) over the attosecond pulse and use 
that no contribution from the first term on the LHS of 
([T]) appears since no continuum electrons are present in- 
side the metal long before or after the pulse. The equa- 
tion describing the transport and energy distribution of 
electrons as a consequence of the attosecond pulse ac- 
cordingly reads 

X(E)(Sl ■ V)N(r, SI, E) + N{r, SI, E) = \{E)P{r, E, SI) 
^ l N(r,Sl',E') i 
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where N(r, SI, E) — dtv<& is the number of elec- 
trons in point r passing unit area along SI with energy 
E, and P(r, SI, E) = J_ dtJ is the number of source 
electrons originating from a single pulse. Equation ([2]) is 
still too complicated to be solved analytical, so to pro- 
ceed we introduce a series of approximations. The mean 
free path X(E) is ~ 5 A in the entire range 30 - 100 
eV [20J, and is taken as a constant A = 5 A in the 
following. The electron distribution within the metal is 
assumed to depend only on the distance to the surface, 
\z\. The skin-depth Sx for the uv radiation is several 
nm's even at grazing incidence, i.e., 5x ^> X(E). Ac- 
cordingly, the spatial derivative in |T]) is also neglected, 
and the source term P is independent of z. The problem 
is then formally identical with the electron beam case 
considered previous, and we follow (l7j and expand N, S 
and P in Legendre polynomials in the angle between the 
velocity of the secondary electron and the surface nor- 
mal. We use that for electron energies up to ~ 100 eV 
the spherical part dominates, i.e., keep only the I = 
component in the Legendre expansion, and we then take 
the s-wave result So = 2/E' [17|, where the factor of 2 
takes into account that electrons may appear at energy 
E either by scattering or by transitions from the con- 
duction band, and where the subscript denotes £ — 0. 
The corresponding Boltzmann equation reads N (E) = 
\P (E) + dE' N {E')-^r, and as is readily seen by in- 
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scrtion, the solution is N (E) = — A J^dE'f^-^ 

We integrate this expression by parts, and find that the 
number of electrons at energy E per attosecond pulse, 
i.e., the photoelectron spectrum, is given by 

/•CO T?l 

N (E) = XP Q (E) + 2 j dE'^XP (E'). (3) 

The first term in ([3]), describes direct electrons and the 
second secondary cascade electrons. 



'-S(E,Sl; E',Sl'), (2) as T, 



To correct for the finite probability of escaping through 
the surface, we multiply Nq(E) from ([3]) by the nor- 
mal component of the free-electron transmission prob- 

ability T{E) = 4^/l-f /(l + yfl - f with V 
the sum of the work function and the Fermi energy, Ep, 
of the metal. 

To evaluate the spectrum, we see from Q that we 
need the number of electrons XPo(E) created with en- 
ergy E, within the distance A of the surface. The ef- 
fects of electron-electron scattering will not be too sen- 
sitive to the exact form of XPo(E), so for convenience 
we follow and consider a T-matrix formal- 

ism for a jellium-like metal. In this model we take 
Xp o = Eoccupicd states | T /»I 2 P(£/X where the T-matrix 
element for the uv-driven transition from the initial 
\^i(t)) and to the final state may be written 



>, - -if^dt^fitWxim^t)) and E f is the 
energy of the final state. The interaction is Vx{t) — 
r ■ eEx(t), with e the polarization vector and Ex(t) the 
field strength of the uv pulse. We include conduction 
and localized core electrons. The latter are modeled by 
&i oc {r) = -^7= J2r e lk ' R tpioc{r — R), where R runs over all 
atom positions within the material, ipiocir) is the atomic 
localized state and N is the total number of atoms in 
the metal. To model the surface, it is required to make 
linear combinations of these solutions, which vanish at 
the surface, z = Q. This is fulfilled by taking the initial 

localized states as ^f oc (r) = 4>f oc {r) — i>i~oc(f) inside the 
metal. Outside the metal, the wave function is vanishing. 

The conduction electrons are modeled with free- 
electron states [lH 14 1 . Inside the metal, ^ b = 

, where V is the total vol- 

F, 



-ik z z 



ume of the metal and 7 = \/2Vq — k' 2 , k 2 /2 < Ep- 

For the final states, we use free-electron states, damped 
within the metal to account for the mean free path, i.e., 



0> f (r,t) 
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"e 1 2 . The density of states 



needed to perform the integration over the continuous 
band of k points within the solid, is with the current 



choice of normalization p(k)d 3 k 
duction and core electrons 22 



T d?k, for both con- 
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When comparing the 
amount of electrons released from the conduction band 
and the core states to experimental results [3j, the con- 
tribution from the delocalized conduction band, which is 
not an accurate description of tungsten[23j, comes out an 
order of magnitude too small. The relative contribution 
has been scaled to represent the experiment. 

We apply our theory to tungsten and take experimen- 
tal parameters Q: We use a gaussian envelope for the 
uv pulse, with 300 asec intensity-FWHM and a central 
frequency 91 eV. We use an experimental Fermi energy 
of 4.5 eV for the conduction band, a hydrogenic localized 
ls-state with a binding energy of 32.5 eV and a tungsten 
work function of 5.5 eV. Note that the series of experi- 
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FIG. 1: (Color online) Calculation of the photoelectron spec- 
trum for a 300 asec uv pulse, including (a) only the direct 
electrons and including (b) both direct and scattered elec- 
trons. The full (black) curve is with the cascade of secondary 
electrons while the dashed (blue) curve includes only one scat- 
tering. The contribution from the conduction band has been 
scaled to be comparable to the contribution from the core 
band. The dotted curve in (b) show data from 



merits [1, 0, & [3] apply an analysis based on laser- induced 
side bands and focus on energies so close to the main uv 
peaks, that we have not been able to assess the influence 
of the background in those data. 

In Fig. [TJ we consider the spectra without assisting 
ir field to clearly display typical spectra of (a) direct 
electrons and (b) secondary electrons that have under- 
gone scattering. Figure QJa) shows that direct emission 
leads to two peaks, representing electrons released from 
the deeper-lying localized states and from the conduction 
band. When the electron-electron scattering is included, 
the two peaks both have a tail of lower energy electrons. 
This situation is shown in Fig. [ljb) where we also show 
the experimental results without ir field [3j . We see that 
the inclusion of a cascade of electron scatterings explains 
the background in the spectrum. As an alternative to 
the cascade theory, we have calculated the electron dis- 
tribution including direct electrons and electrons that has 
suffered exactly one scattering, similarly to 24(. This re- 
produces the amount of secondaries with energy near the 
highest energy direct peak, but fails to account for the 
very many low-energy electrons that are released, sup- 
porting that multiple scatterings and a cascade is at play. 
Note that the experimental data have been subtracted 
for electrons stemming from above threshold ionization 
(ATI) 3] , and hence these do not contribute to the back- 
ground. 

Due to space-charge effects not included in our model, 
the experimental peaks are expected to be upshifted by 
approximately 3 eV [3j. While this shift is seen for the 
core-level electrons in Fig. [ljb), the experimental and 
theoretical peaks pertaining to the conduction electrons 
coincide. This accidental agreement is associated with a 
short-coming of the free electron model which captures 
the qualitative features but does not quantitatively de- 
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FIG. 2: (Color online) Results for tungsten, with a laser inten- 
sity 5 x 10 10 W/cm 2 . (a) Direct emission spectra with scaled 
contribution from the conduction band, (b) Emission spectra 
including scattering to lower energies, (c) Centre-of-energy 
analysis showing the temporal delay of the 4f-electrons (full, 
blue) relative to the conduction band electrons (dashed, red) 
for only the direct electrons, (d) The same as (c) but includ- 
ing scattering. Also shown are experimental data from 
cross: 4f electrons; plus: conduction band electrons. The en- 
ergies of the 4f-electrons has been multiplied by 1.1 in (c) and 
by 1.8 in (d). Notice the difference is scale in (c) and (d). See 
text for laser parameters. 



scribe the tungsten conduction band, which holds a lot 
of structure [23J . Consequently, the center of the peak 
comes out of our calculation a few eV too high, since 
many of the conduction band electrons are actually a 
few eV deeper bound [23j |. 

We now turn to the streaking spectra. We follow [3], 
and invoke the streaking directly in the source term by 
adding the appropriate Volkov phase to the final free- 
electron states. For the ir pulse we use a sine-squared 
envelope, a duration 13 fs, a wavelength of 750 nm, and 
an intensity of 5 x 10 10 W/cm 2 . Comparing Figs. 
and (b) shows the effect of including electron-electron 
scattering as the primary photoelectron produced by the 
attosecond uv pulse propagates to the surface and pro- 
duces a cascade of secondary electrons. We see that the 
structure and magnitude of the background mathes quite 
well that seen in experiment after ATI electrons has been 
subtracted (Fig. 2(b) in H). 

Figures Etc) and (d) show the corresponding centre- 
of-energy (CoE) spectra. The CoE analysis is performed 
on the interval 44 eV - 63 eV (4f-electrons) and above 66 
eV (conduction band electrons). The 4f-electron curve is, 
however, on top of the conduction band electron tail and 
is damped more than the conduction band curve. As in 
[ll ] , the 4f curve has been multiplied by a factor (in our 
case 1.8) to make the 4f and conduction band amplitudes 
comparable. We see from a comparison of Figs. Etc) and 
(d) that the inclusion of the background produced by 
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secondary electrons leads to a reduction of the amplitude 
in the CoE spectrum and a very good agreement with 
the measured spectrum Q. We note that the amplitude 
in the CoE data can be reproduced with the present T- 
matrix theory not including the background, but for an 
unrealistic low intensity of the femtosecond laser pulse. 
In this respect it is essential to have access to data as in 
Fig. [Ib). 

Although the issue of the time-delay is not a prime con- 
cern in the present work, we make the following comment. 
We find a shift of ~ 300 asec which is larger than what 
was found experimentally Q and in a similar model [3] • 
This discrepancy might be due to the failure of the jel- 
lium model [3, EH to describe the structured conduc- 
tion band in tungsten [23[, and points to the need for 
further work to fully understand the origin of the time 
delay, e.g., along the lines of time-dependent wavepacket 
calculations 16J. In this connection, we also note that 
in [3, Q] and this work, the streaking takes place inside 
the metal assuming the penetration depth of the ir field 
to be much larger than the mean free path, while in [l6| , 
the streaking occurs as the electron wave packet leaves 
the surface. 

With the inclusion of electron scattering within the 
metal, we are able to reconstruct the significant tail of 
electrons at lower energies observed experimentally in at- 
tosecond photoelectron spectroscopy from metals Q • In 
monochromatic experiments there is generally relatively 
few electrons scattered to lower energies, as compared to 
what is seen in the attosecond regime. This is due to 
the fact that the amount of electrons scattered is pro- 
portional to the total amount of excited electrons. In 
an experiment with a spectrally broad pulse, the tail of 
electrons becomes comparable to the main peak of direct 
electrons. 

In conclusion, we have presented a simple model for 
electron scattering effects in the attosecond domain for 
electron propagation through a metal, and we have shown 
that the model reproduces two sofar unexplained cen- 
tral features of the spectra in this regime (i) the back- 
ground and (ii) the amplitude in the streaking spec- 
trum. Electron-electron scattering is an important pro- 
cess to consider to fully understand attophysics in the 
condensed phase. The significance of the secondary elec- 
trons relative to the primary electrons is much greater 
in attosecond spectroscopy than in conventional photo- 
electron spectroscopy due to the broad bandwidth of the 
pulse. Our findings ([3|) show that the relative amount 
of secondary electrons near the direct peak scales as 
E" 1 . Hence, to obtain higher energy resolution, it would 
be advantageous to use higher photon energies as may 
become possible with the upcoming free-electron laser 
sources. The present work shows that the classical Boltz- 
mann equation captures important physics in the consid- 
ered regime and, hence, in combination with input from 
quantum mechanical calculations, provides an attractive 



starting point for further theory development. 

After the submission of this manuscript, another work 
using classical transport theory to explain attosecond 
streaking appeared [25(. 

This work was supported by the Danish Research 
Agency (Grant no. 2117-05-0081) 
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